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PATTERN FORMATION (II): THE TURING 
INSTABILITY 

YAN GUO AND HYUNG JU HWANG 



Abstract. We consider the classical Turing instability in a reaction- 
diffusion system as the secend part of our study on pattern forma- 
\ tion. We prove that nonlinear dynamics of a general perturbation 

of the Turing instability is determined by the finite number of lin- 
ear growing modes over a time scale of In j-, where 6 is the strength 
p ^ ' of the initial perturbation. 

< 
■3 

■ 1. Growing modes in a reaction-diffusion system 

In this section we summarize the classical linear Turing instability 
criterion for a reaction-diffusion system. Consider a reaction-diffusion 
>- ■ system of 2-species as 

(1-1) ^ = v-{D,{uy)vu) + f{uy), 

2 : ^ = V ■ (^2 (uy) vv) + g (uy) , 

^ . where U {'x.,t) y {:K,t) are concentration for species, Di,D2 diffusion 

! coefficients, /, g reaction terms. 

^ ! In this paper we consider a d-dimensional box T'^ = (0, vr)'^ , d = 

1, 2, 3, with Neumann boundary conditions for U and V, i.e., 

^ ■ dU dV 

k> ■ (1.2) 7^ = 7^ = at = 0, TT, foTl<i<d. 

rS ' oxi dxi 

Hi - - 

I Homogeneous steady state U = Uy=V forms a steady state provided 

(1.3) = / (uy) = g (uy) . 

In this article, we study the nonlinear evolution of a perturbation 

u{x,t) = U{x,t) - U, v{x,t) = V{x,t) - V 
around [U, V] , which satisfies the equivalent react ion- diffusion system: 

3u — — — — 

(1.4) _ = V ■ (Di (m + U,v + V) Vu) + f {u + U,v + V) , 

Ol) — — — — 

(1.5) _ = V ■ (1^2 (m + U,v + V) Vv) +g{u + U,v + V) . 
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The corresponding linearized system then takes the form 

(1.6) Ut = DiV\ + fuU + f,v, 

(1.7) vt = D2V'^v + guU + g^v, 

where D, = D,{U,V), = D2{U,V), U = ^{U,V)J, = 

We use [■, ■] to denote a column vector, and let 

w{x,t) = [u{x,t),v{x,t)], W = [U,V]. 

Then the original nonlinear system (jl.4|) and (|1.5|) can be written in a 
matrix form: 
dw 

(1.8) — =V-(DVw)+F 

= (DVV+Aw) + ({V ■ {DVw) - DV^w} + F-Aw) 
=C (w) +Ar(w). 

where 

/ (w + W) _ \ Di 

V ^2 (w + W) J ' V ^2 

V (7 (w + W) ; ' \gu gv 
Let q = (gi, .., g^) G ^^ = (N U {0})"' and let 



eq(x) = ]^cos (giOJi) , 

i=l 

where q Gfi. Then {eq(x)}qg^ forms a basis of the space of functions 
in T*^ that satisfy Neumann boundary condition ()1.2j) . 

We look for a normal mode to the linear reaction-diffusion system 
()1.6|1 and ()1.7|) of the following form: 

(1.9) W (X,t) = FqCXp (Aqt) eq(x), 

where rq is a vector depending on q. We substitute (|1.9|) into (|1.6|) - (jl.7|) 

to get 

_ f fu~ Diq^ fv_ 

where = Ylt=i1i- nontrivial normal mode can be obtained by 
setting 

det ( ~ ^^^^ -fv _ 

-9u Aq - + D2q^ 
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This leads to the following dispersion formula for Aq: 
(1.10) 

Xl+{-fu+Diq^-gv+D2q^}K+{{fu - D,q^) {Qv - D2q^)-f,gu} = 0. 

We assume first that without diffusion, the Aq has negative real part 
(stable) : 

(1.11) tr A = /„ + ^„ < 0, det A = f^g^ - Lqu > 0, 

On the other hand, in the presence of diffusion, we assume the following 
diffusion-driven (linear) instability criterion by requiring there exists a 
q such that 

(1.12) (/„ - (^^ _ D2q') _ f^g^ < q, 

which ensures that (jl.lOj) has at least one positive root Aq. 

Remark 1. To satisfy / ti.iij) and M.l^) . the discriminant for the qua- 
dratic equation for q^ in lll.l^} must be positive: 



(1.13) {fuD2 + -g,D,) > 2VD,D2 det A > 



which means the range of inhibition \/ D2I \gv\ is larger than the range 



of activation y-Di/ \fu\- From il.ll]) and M.IS^) . it follows that 

(1.14) fugv < 0, and f^g^, < 0, 

and we have only two cases for A : 

A=(+"^ or A = ^ ^ ^ 



where formal case is called activator-inhibitor (or predator-prey) and 
the latter positive feedback. It also follows from il.ll]) that 

Di ^ D2. 

For given q G f2, we denote the corresponding eigenvalues by A±(q) 
and eigenvectors by r-|-(q). We split into the three cases for the linear 
analysis: 

(1) Generic case where we have two independent real eigenvectors 
and we denote 

f^generic = {q G such that r+(q) 7^ r_(q)}. 

By an elementary computation of the discriminant of ()1.10|) . we have, 
except for only finitely many q, 

{D, - D2) q^-ti A {D, + D2) g'+4 {UD2 + g.D^) g'+(tr Af-A det A > 0, 
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since Di — D2 7^ 0. Therefore, there are two distinct real roots such 
that 

A_(q) < A+(q) 

for large q. Since 7^ in ()1.14|) . the corresponding (linearly indepen- 
dent) eigenvectors r_(q) and r+(q) are given by 



:i.l5) r±(q) 



1, 



fv 
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It is easy to see from ()1.12|) that there exist only finitely many q such 
that A+(q) > 0. We therefore can denote the largest eigenvalue by 
A max > and define 

^max = {q G such that A+(q) = A^ax }• 

We also denote z/ > to be the gap between the Amax and the rest. 
Moreover, there is one (possibly two) having Aq (g^) = Amax when 
we regard A^ as a function of g^. 

(2) Defective case where we have the repeated real eigenvalues and 
eigenvectors: 

Note that there may be possibly one (so finitely many q) such 
that from (irTT|l 

(1.16) A+(q) = A_(q) = A(q) = {/, + 9v - (^i + ^2) q'}/2 < 
and r+(q) = r_(q) = r(q) and we denote 

^defective = {q G such that r+(q) = r_(q)}. 
In this case we find another independent vector 

r'(q) = [0,i] 

Jv 

satisfying {A — A(q)/) r'(q) = r(q). 

(3) Complex case where we have complex eigenvalues for q and we 
denote it by i^compicx = ^ — (^^gcncric U f^defective)- For q G ^^compicx, we 
denote A+(q) = ReA(q) + ilmA(q) and r+(q) = Rer(q) + ilmr(q). 
Then we have A_(q) = ReA(q) — ilmA(q) and r_(q) = Rer(q) — 
i Imr(q). Notice that Re A(q) < as in ()1.16|) . and Rer(q) and Imr(q) 
are linearly independent vectors. 
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Given any initial perturbation w (x, 0), we can expand it as 
w (x, 0) = ^ Wqeq(x) = {Wqr-(q) + w+r+(q)}eq(x) 

qGf^ q^^gcncric 

+ 5Z {^qr(q) +^^qr (q)}eq(x) 

q ^ ^defective 

+ Yl {<'Rer(q) + ^S-lmr(q)}eq(a;), 

q ^ ^complex 

so that 

(1.17) Wq = w~r_(q) + w+r+(q) for q G ^generic, 

Wq = Wqr(q) + w'^v'{q) for q G ^^defcctive, 
Wq = wj'' Re r(q) + tyj^ Im r(q) for q G f^compicx- 

The unique solution w (x, t) = [u (x, t) , t> (x, t)] to ()1.6|) - ()1.7p is given 
by 

(1.18) 

w {x,t) = Y exp (Aqt) + w+r+(q) exp (A+t)}eq(x) 

q Gfigcncric 

+ 5Z { (^^qi'lq) + ^qr'(q)) + w'^r{q)t} exp (Aqt) eq(x) 

efective 

+ Y {Wq^(Rer(q)cos[(lmAq)t] -lmr(q)sin[(lmAq)t]) 

q Gf^complcx 

+ (Re r(q) sin [(Im Aq) t] + Im r(q) cos [(Im Aq) t])} exp[(Re Aq) t]eq(x) 
= e^*w (x, 0) . 

For any u (•,t) G [L2(T"')]^ we denote ||u(-,t)|| = \\u{-,t)\\^,. Our 
main result of this section is 

Lemma 1. Assume that il.ll]) and the instability criterion il.l^) are 
valid. Suppose 

w (x, t) = [u (x,t) , V (x,t)] = e^^w (x, 0) 

as in m.l^} is a solution to the linearized reaction- diffusion system 
/ li.6|) - P^ with initial condition w (x, 0). Then there exists a constant 
Ci > 1 depending on U, V, Di, D2, A such that 

||w(-,t)|| < Ciexp(A™axt) ||w(-,0)||, 

for allt > 0. 
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Proof. We first notice that from tlie quadratic formula for , for q 

large, 

|det[r„(q),r+(q)]| = ^^Z^ > J^'K '- 

\jv\ 

Thus solving ()1.17|) yields, due to Di 7^ -D2, 
'I det[r_(q),r+(q)] 

<C|Wq|, 

Since Aq < 0, for q G ^^defective , we have 

texp(Aqt) < C. 

Moreover, recall ReA(q) < for g G ^^compicx- Thus we deduce the 
Lemma on the linear growth rate by the formula ()1.18|) . □ 



2. Main Result 

Let 6* be a small fixed constant, and Amax be the dominant eigenvalue 
which is the maximal growth rate. We also denote the gap between the 
largest growth rate A^ax and the rest by z/ > 0. Then for 5 > arbitrary 
small, we define the escape time by 

(2.1) ^ = 5exp(A^,,T^), 

Theorem 1. Assume il.ll}) and that there exists q^ = 'Y2'i=iQi satis- 
fying instability criterion il.l^) . Let 

Wo(a;) = ^{u;qr_(q) + u;+r+(q)}eq(x) 

+ Yl {^qi"(q) + ^qi^'(q)}eq(a;) 

ci^^ defective 

+ Yl {w^J'Rer(q) +w5"Imr(q)}eq(x). 

ci^^ complex 



or equivalently 



Our main theorem is 
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G such that ||wo|| = 1. Assume Di,D2 f\g E C"^ near W, so that 

there exists rj > 

(2.2) 

2 

C,= max {y^\\D,{W+w)\\c2 + \\f{W+w)\\c2 + MW+w)\y<^. 

Then there exist constants 6o > 0, C > 0, and 9 > 0, depending on 
U, V, Di, D2, f, g, such that for all < 6 < 60, if the initial perturbation 
of the steady state [tJ, V] in il.!^) is 

(x,0) = (5wo, 
then its nonlinear evolution w^{t,x) satisfies 

(2.3) \\w\t,x) - 66^"'^^^ w+r+(q)eq(x)|| 

max 

for < t < T^, and u > is the gap between Xmax Ojf^d the rest of Re Aq 

m (ma). 

We notice that for < t < T^, 5e^'"="'* < 9, is sufficiently small. The 
initial proffie wq is any function. In particular, as long as u'q^ 7^ 
for at least one qoGf^max (generic for a general perturbation), the 
part of its fastest growing modes satisfies 

||^gA„,a.t J2 w+r+(q)eq|| > (5e^-^''*|w+ ||r+(qo)|, 

which has the dominant leading order of 5e^™="'*. Our estimate ()2.3|) 
implies that the dynamics of a general perturbation can be character- 
ized by such linear dynamics over a long time period of eT^ < t < T^, 
for any fixed constant e > 0. In particular, choose a fixed qo G fimax 
and let 

wo{x) = ] — r^^io a; 



then if t = T'', 



hence 

||w''(t, ■)\\>9- Cj^'^/^-- + > 9/2 > 0, 

which implies nonlinear instability as 5 — >^ 0. The instability occurs 
before the possible blow-up time. 
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Reaction-diffusion systems are often employed to study chemical and 
biological pattern formation and have received much attention from sci- 
entists jS], [13], jlSj, [H]; since the pioneering work of Turing jTT] 
in 1951. This symmetry breaking instability is called diffusion- driven 
instability, since the presence of diffusion and the difference of diffusion 
coefficients are essential for the instability mechanism and nonuniform 
pattern formation. After some experimental results such as in |2], 
|15j . more extensive and serious works began towards this Turing-like 
pattern formation across many fields of study. Our result can be inter- 
preted as a mathematical description of early pattern formation. Each 
initial perturbation can be drastically different from another, which 
gives rise to the richness of the pattern; on the other hand, the finite 
number maximal growing modes determine the common characteris- 
tics of the pattern, over the time scale of ln|. In comparision with 
an earlier different result along this direction ^Hj: First of all, the 
reaction-diffusion system considered here is not scaled. Secondly, our 
initial perturbation is more general, need not be close to the space of 
finite number of maximal growing modes. Thirdly, a precise estimate 
of the time scale (In |) for pattern formation is given here, without 
an a-priori assumption for the smallness of the perturbation later in 
time as in jTH|- Lastly, based on Guo-Strauss' bootstrap argument, our 
proof is much simpler and direct. 

3. Bootstrap Lemma 
We state existence of local-in-time solutions for (jl.4|) - (|1.5|) . 

Lemma 2. (Local existence) For s > 1 {d = 1) and s > 2 {d = 2,3), 
there exist aT > and a constant C depending on U ,V , Di, D2, f, g 
such that ||w(t)||j^s is continuous in [0,T), and 

||w(t)||^. <C||w(0)||^.. 

We now derive the following energy estimates for ci-dimensional reaction- 
diffusion system with d = 1,2, 3. 

Lemma 3. Suppose that [u {x,t) , v {x, t)] is a solution to the full system 
ra-Ca;. Then for \\w{t)\\H2 <v, 



1 d 
2dt 



V [ {\duf + \dvf}dx 

+ E / [^\^du\' + D,\WdvAd^+^Yl I \M 
< CoCi||w||h2||V^w|P + C2|b||^ 
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where Co is the universal constant while Ci = CoCr,{l + rf) and 

,2 \ 3 



Co 



2\gv\ 



m 

Proof. We first notice that the reaction-diffusion system p.4j) - ()1.5|) pre- 
serves the evenness of the solution w{x,t), i.e., if w(x, t) is a solution, 
then w(— Xj, t) is also a solution. We can regard the Neumann problem 
as a special case with evenness of the periodic problem by standard 
way of even extension w{x,t) with respect to one of the Xj. For this 
reason we may assume periodicity at the boundary of the extended 
periodic box 2T^ = (— vr, vr)'^. Since now there is no contributions from 
the boundaries, we can take second order 9-derivative of (jl.8|) to get 

(3.1) ~[ \dw\^=[ dw'^dC{w)+ [ dw'^d^iw). 

We first treat the last nonlinear term: 

- / {VdwY[d{D (w + W) Vw} + DVdw] + {Vdw^d (F-Aw) 

<C||d(w + W) - D\\^\\Vdwf + C\\{VD) (w + W)||^||Vw||^||9w|| ||V9w|| 
+ C\\{dD) (w + W)||^ ||Vw||^4 ||Vw||^||vaw|| 

+ C\\{d¥) (w + W)||^ ||Vw||^ ||Vw|| ||V9w|| +C||VF (w + W) -A\\^ \\dw\\ ||V5w|| . 
We apply the following the Sobolev imbedding to control ||w||oo 

(3-2) llfi'llL°°(2Td) ^ CQ\\g\\^2(^2T<i) ^ 

for (i < 3. Moreover, from the periodic boundary conditions. 



Vu= i = 0, 

2T'' J2T'* 

we also use the Poincare inequality 

(3.3) \\g\\<\\g\\Lii2T'i)<CQ\\Vg\\ ifrf<3. 

to further get 

||Vw|U < Co ||Vw||^, <CoY. PVwII. 

\d\=2 

where Cq is a universal constant. From (j2.2j) and the assumption 
||w||iy2 < r/, the last nonlinear term in ()3.H1 is bounded by 

CoC^(l + r/)||w||^2 ||V9wf . 
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We now estimate the second quadratic term in (jH.lj) 

-[ {DilVduf + D2\Vdvf} + g, [ \dv\ 

+ {fv + gu) / dudv + fu j \du\^ . 
The last two terms are bounded by 

[fv + 9u) / dudv + fu \duf 



,2 



J2T' \9v\ J2T^ 

Thus we can bound the hnear term in ()3.H) by {jjv 

-[ {Di\Vduf + D2\Vdvf}-^-^[ \dvf 

+ {^h^ + fu}f \du\\ 

^ \9v\ J2T<^ 

By the interpolation between ||V5m|| and , the last term 
above is bounded by 

for any a > 0. We can choose a such that 



^ 2 1^,1 '"^ 2 
Collecting terms, we conclude the proof. □ 

We are now ready to establish the bootstrap lemma, which controls 
the growth of w(x,t) in term of its growth nonlinearly. 

Lemma 4. Suppose that w{x,t) is a solution to the full system 
fOl) such that forO<t<T 

\\w{-,t)\\H2 < mm <^ T], 



2C(jCi CqCi 
and 

(3.4) l|w(-,t)|| < 2Cie^-^''*||w(-,0)||, 

then we have for <t <T 

l|wW||^.<C73{||w(0)||^. + e2^-*||w(-,0)|p} 
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where C3 = C? max|^, 1} > 1. 

Proof. It suffices to only consider the second-order derivatives of w(a;, t). 
From the previous lemma and our assumption for ||w||/^2, we deduce 
that 



, , I |2 
2\\U\\ ■ 



\a\=2 ' 

So that by ()H.4|1 and an integration from to t < T, we have 



/ {\du{t)\' + \dv{t)\'] 

r An n2 

\ {|a«(0)p+|9.;(0)p} + 



An nl 

< > / Wdum^ \\'dv[^)?^ \ / ^ e^^— *||w( 

An 



>max 



=2 • 

Thus our lemma follows. □ 

4. Nonlinear instability and pattern formation 

We now prove our main Theorem 1: 

Proof. Let w'^ (x, t) be the family of solutions to the reaction-diffusion 
system (jl.4j) - (|1.5p with initial data w*^ (x, 0) = Swq. Define T* by 

T* = sup |t I ||w^(t) - 5e^*wo|| < Y^exp (K^J)^ . 

Note that T* is well defined. We also define 

T** = sup It I ||w(t)||//2 < min li] ^ ^ 



2CqCi CqCi 

We now derive estimates for norm of w''(x, t) for < t < 
min{r*,T**}. First of all, by the definition of T*, for t < T* and 
Lemma IT] 

3C 

11^(^)11 < exp (Amaxi) • 

Moreover, using Lemma 4 and applying a bootstrap argument yields 

(4.1) ||w'(t)||^, < v^{(5||wo|U2 + 5e^— 

We now estimate the norm of w^{x,t) for < t < min{T*,T**}. 
We apply Duhamel's principle to obtain 

(t) = (5e^*Wo - / e^(*-")7V (w^ (r)) dr, 
Jo 
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Using Lemma 1, (jS21), (jSISI), and Lemma 4 yields, for < t < min{T*, T**} 
\\w^ (t) -5e^*wo|| 

<Ci [ e^"^^^(*"^) ||{V ■ (DVw^) - DV'^w^} + F-Aw^\\ dr 
Jo 

Jo 

+ Ci [ e^-^^^^^-'^llZ^IIci II Vw^ 
Jo 



:^)|L4||Vw^ i^niL^dr 



+ Ci [ e^--(*~")||F||c2 ||w'^(r)||^ ||w'^(r)||rfr 
Jo °° 

<CiCo'C, Te^— ||w^(r)||^,rfr. 
Jo 

from assumption ()2.2p with ||w||iy2 < rj. We plug ()4.ip with t = r to 
further obtain 



(4.2) ||w''(t) -5e"wo|| 

Jo 



We now choose 6' in in ()2.ip to satisfy 
(4.3) C',CsC,e < 



4 ' 

5i D2 



(4.4) < min <^ r/, 



2C0C1 CqCi 

We now prove by contradiction that for 6 sufficiently small, 

< min{T*,T**}, 



and therefore our theorem follows from ()4.2|1 . by further separating 
q £ ^max and move q ^ l^max in to the right hand side . 

If T** is the smallest among T\ T* and T**, we can let t = T** < 
in (jHH) 

||w^(T**)||^, < v^{5||wo||^2 + 5e"'} 

= ^{<5||wo||^^2 + ^^}<2v^^. 
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for small 5 such that (5| |wo| li^a < 9. By the choice of 9 in ()4.4j) . we have 
||w(T-)||^.<min{r,,^,^}. 

This is a contradiction to the definition of T**. 

On the other hand, if T* is the smallest among among T\ T* and 
T**, we can let t = T* in (jOI) to get 

llw'' (T*) -(5e^*wo|| 

1 



^max 

|2 



'^max ^max 

C\ r- \ T* 

< —6e , 

for CgCsC j'^l'f''^ < 1/4 for 6 small, by our choice of 9 in (jOD- This 
again contradicts the definition of T and our theorem follows. □ 
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